A propensity criterion for networking 
in an array of coupled chaotic systems 
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We examine the mutual synchronization of a one dimensional chain of chaotic identical objects 
in the presence of a stimulus applied to the first site. We first describe the characteristics of the 
local elements, and then the process whereby a global nontrivial behaviour emerges. A propensity 
criterion for networking is introduced, consisting in the coexistence within the attractor of a localized 
chaotic region, which displays high sensitivity to external stimuli, and an island of stability, which 
provides a reliable coupling signal to the neighbors in the chain. Based on this criterion we compare 
homoclinic chaos, recently explored in lasers and conjectured to be typical of a single neuron, with 
Lorenz chaos. 

PACS numbers: PACS: 05.40.-a, 05.45.-a, 05.45.Xt 
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An open problem in science is how to build a semantic 
network on minimal assumptions. In a living brain an 
external stimulus localized at some input spreads over a 
large assembly of coupled neurons building up a collec- 
tive state univocally corresponding to the stimulus. The 
current conjectures and the preliminary experimental ev- 
idence 0, on this time dependent networking problem 
lead to a new paradigm wherein perceptions require mu- 
tual synchronization of neuronal spikes. A model of dy- 
namical encoding by networks of competing neurons has 
been recently introduced Q, however the core issue of 
synchronizing a large number of neurons as it appears 
from experiments [2| has not been addressed thus far. 

With reference to this issue we introduce a minimal 
model built upon simplicity requirements, namely we 
consider linear and symmetric interneuron coupling; fur- 
thermore, we take just nearest neighbor coupling, avoid- 
ing architecturally complicated connections. Under these 
assumptions we address the question of how complexity 
arises, calling complexity the fact that a composite sys- 
tem displays collective properties not directly deducible 
from the dynamical behavior of its constituent elements 
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Here we focus on the collective response of an array of 
coupled dynamical objects to a localized external stim- 
ulus, and provide a propensity criterion for networking, 
that is, for organizing in a collective state univocally re- 
lated to the stimulus. Having in mind time dependent 
networking, as it is the case of biological communication, 
we examine mutual synchronization of a one dimensional 
chain of identical objects in the presence of a stimulus 
applied to the first site. We first describe the character- 
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istics of the local elements, and then the process whereby 
a global nontrivial behaviour emerges. 

Two independent requirements must be fulfilled by 
each element. First, it should have a sensitivity region in 
order to easily respond to the neighbour coupling; second 
it has to provide a strong enough signal to relay the in- 
put stimulus along the whole array. The first requirement 
suggests to recur to dissipative chaotic systems; in fact 
modifying a regular individual dynamics confined to a 
stable attractor would imply a consistent expenditure of 
energy and time [f| , whereas sweeping through the man- 
ifold of unstable periodic orbits which make up a chaotic 
attractor is a fast cost-less operation 0. The second re- 
quirement is conflicting with the previous one; in fact, it 
implies islands of stability within the chaotic orbit, out of 
which to extract a reliable driver for the next neighbour. 

The twofold problem is solved recurring to large spikes 
emerging out of a small chaotic background; indeed a 
weak inter-site coupling will provide a discrete synchro- 
nization associated only with the large spikes, the chaotic 
background being not effective. This discrimination will 
amount to the approximation inducing the transition 
from an individual to a collective description. At vari- 
ance with the standard chaotic synchronization scheme 
where two identical systems synchronize along the 
whole orbital evolution, here the synchronization occurs 
because a large spike of one system is forcing the neigh- 
bour to escape away from its chaotic region, thus yielding 
its own spike 0, . 

Among the chaotic systems, those best suited to the 
emergence of a new hierarchical level should thus be char- 
acterized by temporal windows of stability and chaos 
within each orbit. 

We give substance to these considerations with refer- 
ence to the heteroclinic transfer forth and back between 
a saddle focus (SF) and a saddle node (SN) under the so 
called Shilnikov condition Such a behavior has been 
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explored both experimentally [ll| an d theoretically 0] 
with reference to a CO2 laser with feedback. For sake 
of brevity it will be called HC (homoclinic chaos) [l3j| . 
In fact the mere homoclinic return to SF would provide 
a chaotic transient 0] but not assure a regolar motion 
away from SF ; on the contrary, the further presence of 
SN yields the stability island necessary for networking. 

The dynamics is characterized by a sequence of spikes 
with widely fluctuating time intervals T. Such a struc- 
ture underlies spiking behavior in many neuron [T5I , 
chemical laser [u|, and El Nino [3] systems. It is 
important to note that this dynamics is highly nonuni- 
form, in the sense that the sensitivity to small perturba- 
tions is high only in the vicinity of S along the unstable 
directions. A weak noise thus may influence T signifi- 
cantly 0. . 

In the system we consider |l2j the chaotic behavior is 
confined in a small neighborhood of SF where it fulfills 
the Shilnikov condition 7 < a ,7 being the real part 
of the expansion rate on the unstable manifold of SF and 
a the contraction rate on the stable manifold of SF. Pre- 
cisely the linearized dynamics around SF is ruled by the 
leading eigenvalues (— a,7±ia->) with a = 7.58, 7 = 3.52, 
all the other eigenvalues having very large negative real 
parts. The exiting trajectories along the unstable man- 
ifold reenter the stable one after a large orbit in phase 
space, corresponding to the heteroclinic approach to SN. 
The phase space orbit appears as a confined chaotic tan- 
gle, concentrated around the saddle focus which is con- 
nected to a wide regular section, the two parts closing 
the orbit (Fig. 1-left). 

For such a system synchronization has been demon- 
strated against an external clock [l^, or against a pre- 
vious time slot of itself, presented after a suitable delay 
[2oT| or after a low pass filter |2l|. Furthermore, HC is 
robust against noise 0|. As we go from one system to 
an array of coupled identical systems 0, E| , mutual syn- 
chronization occurs either spontaneously or as a response 
to an external forcing applied to a single site. 

We recall that the coup ling is realized by replacing 
in one of the HC equations 12] a scalar component x\ 



(i = site index) with x\ [1 + e(x\ 



(+1 



2 < x\ >)] 



where < x\ > denotes a moving time average. The range 
of coupling strengths e considered here is between and 
0.25. 

To show the advantage of HC in comparison to more 
conventional types of chaotic behaviour, we consider the 
well known Lorenz model In terms of fixed points, 
also the rectified Lorenz system(i.e.,the Lorenz model 
plus an inversion operation around the origin) is charac- 
terized by one SF and one SN; however for the standard 
values of the control parameters b=8/3,<r=10, r=28, the 
eigenvalues around SF are (—a ± iuj, 7) with a=14.48 
and 7 = 0.4119, thus very far away from the Shilnikov 
condition. 

In fact, morphologically the two cases appear very dif- 
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FIG. 1: Comparison between HC(homoclinic chaos) [Tll| (left 
column) and Lorenz chaos (right column). The top row is a 
phase space projection over two dynamical variables; SF de- 
notes the saddle focus and SN the saddle node, on the left the 
two SF map one onto the other after an inversion (a? — > —x, 
y —* —y) around the origin. The intermediate row shows the 
time series for variables x\ of HC (it represents the laser in- 
tensity in the case of the CO2 laser) and x of Lorenz; in the 
former case, a suitable threshold cuts off the chaotic back- 
ground; in the latter case, no convenient region for threshold- 
ing can be isolated. Bottom row, left : after threshold, the 
new variable S(t) alternates spikes with flat regions where the 
systems has a high sensitivity and short refractive windows 
where the intensity S(t) goes to zero. 



ferent; thus, even though a homeomorphism should map 
one over the other, introduction of extra operations, such 
as a sensitivity threshold, renders such a transfer impos- 
sible, as one can see by comparing left and right columns 
of Fig. 1 before and after thresholding (intermediate and 
bottom rows). 

Precisely, the time plot of one of the HC variables con- 
sists of a train of identical spikes, separated by a vari- 
able interspike interval (ISI) filled with a rather small 
chaotic signal confined within a stripe thinner than 1/5 
of the spike height. As we assemble a large number of 
such dynamical systems in a network and identify an op- 
timal mutual coupling, which yields an efficient transfer 
of the input information along the network, a reduction 
by a factor 3 of that coupling makes already the trans- 
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FIG. 2: Coherence parameter R\% for a driving signal con- 
sisting of a 1 % periodic perturbation of a control parameter 
normalized to Rf ree ( the R value in the absence of pertur- 
bation), plotted versus the normalized distance of the pertur- 
bation frequency from the natural value u>o. Circles : HC; 
squares : Lorenz. For HC the ratio is about 30 at loo 

and it increases up to 10 4 for u twice or more the natural 
frequency; indeed, since synchronization means forcing HC 
away from the SF region, frequencies higher than the natural 
one are easier to synchronize, whereas for smaller frequencies 
HC may have a spontaneous escape from SF. To smooth the 
plot, a small amount of noise (0.5 %) has been added. For 
the Lorenz case, the indicator is always at 1. 



fer inefficient, as shown later in Fig. 3; thus a natural 
thresholding is provided by the alternation large-small 
which makes the bottom plot on Fig. 1 left equivalent to 
the intermediate one for networking purposes, without 
need to apply the threshold as an artifact. 

The emerging spikes display a chaotic time occur- 
rence, which is epitomized by the correlation properties 
of I Sis. Mathematically, this thresholded chaos is de- 
scribed as S(t) — S J2i<i{t — T i) where t; are the chaotic 
spike positions and (ISI)i = tj — ti-\- In presence of 
an external driving signal, the spikes can synchronize 
to it |13| : in case of many coupled systems of this type 
located on an array, wide parameter ranges are found 
within which the individual sites mutually synchronize 
their own spikes, so that the space-time plot of the spike 
positions at each site appears as a regular fabric d, • 

Such an easy mutual synchronization as the response 
to a localized input represents a semantic property ab- 
sent in the Lorenz case where there is no apparent 
scale separation where to consider a threshold. We thus 
identify the propensity criterion with the presence within 
each orbit of two very different amplitude scales, that of 
the large spikes and that of the small chaotic background, 



FIG. 3: Propagation along a chain of N = 20 sites of the 
coherence value Ri of the first site exposed to an external 
stimulus, a) For a fixed input frequency close to the natural 
one coo and different coupling strengths e. b) for e — 0.2 and 
different frequency ratios u/u)q. 



the mutual coupling between sites being operated by the 
spikes whereas the background represents the high sensi- 
tivity region within which the coupling takes place. 

As indicator of successful networking we take the co- 
herence parameter 



R 



< ISI > 
6ISI 



(1) 



where SISl is the square root of the ISI variance (51 SI = 
'< (ISI— < ISI >) 2 >). R has been used to measure 
the amount of synchronization of a single system to a pe- 
riodic stimulus; when synchronization propagates along 
a chain, R can be measured anywhere since it has almost 
the same value on all sites. 

Notice that standard chaotic synchronization can prop- 
agate along a chain of generic chaotic systems, however, 
in general there is no propensity of the first site to syn- 
chronize to a weak external input, as shown by the cor- 
responding indicator R (Fig. 2) for the systems of Fig. 1, 
either in the case of propensity (left) and no-propensity. 
The stimulus consists of a sinusoidal perturbation of fre- 
quency u> applied to a control parameter in the first site 
of the chain, with an amplitude of 1% of the unperturbed 
parameter value. The R value is about 30 for to at the 
natural frequency of the system — 2tt / < ISI >, and 
increases up to 10 4 for larger frequencies, whereas for the 
Lorenz case it is consistentely R = 1 (no coherence). We 
now address the crucial question how the coherence R\ 
induced on the first site i = 1 propagates along the array, 
for different coupling strengths e and frequencies to of the 
input signal. As shown in Fig. 3, reducing the coupling 
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from e = 0.25 to e = 0.08 reduces R by three orders of 
magnitude, thus showing the natural thresholding effect 
occurring in a network of HC systems, without having 
to explicitely take care for the operation leading from 
Fig. 1 (left-middle) to Fig. 1 (left-bottom). 

In conclusion, we have introduced the notion of se- 
mantic network as an array of coupled identical chaotic 
systems which assume a collective state in presence of 
a localized periodic stimulus. The propensity criterion 
appears morphologically as a confinement of the chaotic 
tangle within a small region of the total attractor; the 
corresponding indicator is a tremendous increase in the 
R parameter near the natural frequency. 

At variance with complete synchronization Q , here we 
consider only synchronization of large spikes, intercalated 
by intervals not effective in acting over the neighbours, yet 
highly sensitive to external signals. This split of the dy- 
namical orbit into two different regions is the condition 
to build a nontrivial collective state as response to a lo- 
calized stimulus. 

Two further relevants issues will be dealt with else- 
where, namely i)how the network organizes in presence 
of two stimuli of different period applied to opposite sites; 
and ii)learning of a complex pattern, coded by a non pe- 
riodic sequence of spikes limited in time. 
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